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Abstract 

In a macroscopic realm, in which photons are too many for being counted by any photon counting detec- 
tor, photon statistics can be measured by using detectors simply endowed with linear response. We insert 
one of such detectors in a conventional photon-counting apparatus, which returns a voltage every time the 
detector responds to light by generating a number of elementary charges via its primary photo-detection pro- 
cess. We only assume that, when a single charge is photo-generated, the probability density of the voltages 
is a distribution that is narrow with respect to its mean value. Under this hypothesis the output voltages can 
be suitably binned so that their probability distribution is the same as that of the photo-generated charges, 
that is, of the detected photons. 
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I. INTRODUCTION 



Measuring photon statistics is a useful approach to understand the behavior of any system that 
includes electromagnetic radiation as a part. The investigation of such systems may pertain to 
physics, from astronomy to physics of the matter UJ, as well as to other natural sciences, for in- 
stance biology M\- The availability of photon-counting detectors and methods suitable for any 
situation as to spectral and intensity characteristics of the light to be measured would then be ex- 
tremely desirable. The coverage of the most different spectral ranges is a goal that is pursued by 
the search of novel primary photo-detection processes, including thermal processes occurring at 
cryogenic temperatures. Among the detectors that operate, in essence, as microcalorimeters we 
mention a superconducting transition-edge sensor (TES) with tungsten as the active device ma- 
terial that was recently demonstrated to work as a photon-counter endowed with almost unitary 
quantum efficiency from UV-vis to telecom wavelengths yD. However, we recognize that detec- 
tors based on quantum interactions between photons and sensitive material are largely more used 
than thermal detectors for measuring photon statistics. Detectors based on either external primary 
processes (e.g. electron photo-emission) or internal primary processes (e.g. photo-generation of 
carriers by either photo-voltaic or photo-conductive effects) ensure reasonable values of the de- 
tection quantum efficiency, r/q in the visible and near-IR spectral ranges. The main difficulty that 
still remains with these detectors is that of measuring photon statistics when the charges photo- 
generated in the samples are too many to be counted. Among photoemissive detectors only few 
produce distinct outputs when the number of photo-electrons, m, changes by a unit. The best ones 
are photomultiplier tubes (PMT's) [|5i |6|] and hybrid photo-detectors (HPD's) [16|] that can count 
up to m ^ 5. Relatively more numerous are the photo-emissive detectors that are endowed with 
sufficiently high and sharp gain to provide a sizeable charge in the anodic pulse output for m > 1 
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definitely distinguislied from that for m = 0. Tliey are PMT's available since the 50's that were 
used for the first measurements of light statistics [7, 8, 9, 1^. Nowadays single-photon detectors 
exist that are based on the most different primary photo-detection processes and offer a remedy to 
the lack of good photon counters. In fact, the light to be measured can be split either in space or in 
time prior to detection so that at most one photon at a time hits the detector sensitive area. However 
it must be recognized that these techniques invented for counting photons with intensified CCD 



cameras [11] and multi-pixel and/or position sensitive single-photon detectors IL 
(spatial splitting) or single-photon avalanche photodiodes (temporal splitting) [[16 



14, 



1711 are rather 



cumbersome. Their adoption is only justified by the impossibility of performing direct measure- 
ments with photon counters when the number of detected photons becomes macroscopic. 

The work described here concerns the direct measurement of the detected-photon statistical 
distribution, P^, and is motivated by the fact that, in many of the systems for which measuring 
photon statistics is relevant, artificially lowering m is not permitted either by attenuating the light 
or by shortening the measuring time, Tm- This is the case of fields that modify their properties 
upon attenuation and, obviously, of pulse fields in which Tm cannot be shorter than the light pulse 
duration. It is worth noting that measuring photon statistics when m > 1 in Tm is a problem that 
has been faced since the 60's. In particular Arecchi et al. [ISO suggested a "linear method" in 
which the PMT anodic charge corresponding to the photons detected in Tm was recorded. More- 
over these authors demonstrated that calculating the moments of the statistical distribution of this 
charge and those of the single-electron response (SER) distribution allows obtaining the moments 
of Pjn- Such a result has been used to verify the agreement with the theoretical Pm moments up 
to second order [19]. However, using it to recover Pm would be at least cumbersome owing to 
the need of accurate evaluations of SER distribution moments at any order. We will show that any 
detector based on either an internal or an external primary photo-detection process and endowed 



with two properties rather commonly encountered, allow measuring Pm in a macroscopic realm 
in which photon-counters do not exist. The two properties are: (i) the detector response must be 
linear up to the maximum m of the measurement; (ii) the response for m = 1 must produce a 
standard deviation of the output data that is sufficiently smaller than the mean value. We further 
specify that the detector can be endowed with an internal gain. 



II. MODEL 



With the help of Fig. [T] we first examine how the detector output is processed in a typical direct 
statistical measurement. Normally it is amplified and integrated over a temporal gate, whose dura- 
tion fixes the value of the measure time Tm when a continuous wave light is to be measured. In the 
case of pulsed light, the gate is synchronous and covers the Tm interval in which the current output 
pulse of the detector occurs. The signal is sampled and digitized afterward and converted to a 
voltage V. As indicated in the figure, we represent the overall m-to-v conversion by a single factor, 
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FIG. 1 : (Color online) Measuring apparatus. 



7. Here we will explicitly take into account the statistical distribution of the probability density 
of this conversion factor. The left-hand side of Fig. [T] illustrates the link between the probability 
density of detecting m photons, Pm, and that, P„, of having n photons in the field. Considering the 
effects of the optics that delivers the light to the detector and of the detector quantum efficiency 
r^q < 1 and representing these concomitant losses by an overall photon-detection efficiency, rj. 
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(1) 



n=m 




Obviously in any experiment the value of r] is up-limited by the product of rjq times the coupling 
efficiency of the optical delivery system, but can be diminished at will if filters are inserted into 
the system that delivers the light to the detector. We point out that, as we deal with direct statistical 
measurements, neither the delivery optics includes a fiber looping beam splitter nor the detector is 
a position sensitive one. 

Our aim is to recover Pm for an arbitrary P„ starting from the only experimental data available, 
which are the v voltage values recorded for an ensemble of measurements performed with given 
1] by using an apparatus characterized by a conversion factor 7 with probability distribution py, 
mean value 7 and variance a^. In the following we indicate by Py the probability density of 
the V variable. We assess that we can "measure" P^ as the distribution that we would obtain by 
casting the experimental v values of an ensemble of measurements into a histogram normalized 
to its integral. For ease of writing we represent the bin width by dv, though the variable v is our 
digitized output. The zero of the v scale of the "measured" P„ is set to be equal to the mean value 
of the distribution recorded in a separate experiment performed in the absence of light. As the 
events of having different values of detected photons (i.e. elementary charges generated by the 
primary photo detection process) are mutually exclusive we can write 



in which Pl,*'^ is defined as the probability density of the voltage values v^'^^ that are recorded in 
the events with k detected photons. In Eq. F^^ is the probability distribution measured in 
the absence of light, for which we remind that J^fPl^^df = 0. We note that P^ obviously 
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^ m=0'- v 771=1-'- 1} m=2-'- u • • • 




(2) 



k=0 
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reproduces Pm if the measuring apparatus has photon-counting capability flSL |60 whereas, when 
the Pi'^'*'s do not lead to separate peaks in P^, Eq. Q seems to be useless to reconstruct Pm- The 
latter is exactly the case examined in this paper. 

We consider the central moments Uriv) = {(v — {v)Y) corresponding to the experimental P^ 
and try to relate them to the jJir{fn) = {{m— {m)y') central moments corresponding to the unknown 
Pm probability density. By using properties of the P^'^^ distributions to express /ir(f ), we will find 
relations to /Ur(m) that provide a method to reconstruct Pm- As fI^"* can be identified with the 
probability distribution of the conversion factor 7, we obviously have v^^^ = 7. Owing to the 
hypothesis that the detector response is linear, detecting k > 1 elementary charges corresponds to 
the occurrence of independent events, thus v^'^^ = Yli=i 7«' which all are distributed according 
to p^, and Fi''^ = F^v^ * F^^ * ... * F^^ for k times. Thus we can exploit the following property 
of the cumulants 

k 



E 4"^^ • (3) 

i=l 

By reminding that the lowest order cumulants are = {x), k'^^ = fi2{x), k^^'' = jJ^six), 
K-^^ = /i4(x) — 3[/i2(x)]^, and k^^^ = ^^{x) — 10/i2(x)/i3(x), for the cumulants of the con- 
version/amplification factor we find: h!'^'^ = 7, ^2^'^ = a^, k^^^^ = jl^, k,^^'^ = jl^ — So"^, 
Kg^''' = jl^ — 10(7^/i3, being /ij, the values assumed for the central moments of p^. 

We start by using Eq. ([3]) with r = 1 and then Eq. Q to calculate the mean value of v: 

+00 

{v)=^Y.^Pm=k = {m)^. (4) 

A;=0 

We now calculate the /ir-(f ) moments by applying Eq. 

Y.Pm=k / {v-{v)rF^^^dv 



+ 00 



k=0 

+00 



J2Pm=kMv^'^), (5) 



where we used fir{v^^^) = J^^{v — {v)YFi^^dv. 

For r = 1 Eq. ([5]) obviously vanishes. For r > 2 we make use of the binomial expansion 

/ \ 



w = E 

j=0 



(6) 



VJ 



which, once substituted into Eq. ([5]) and using Eq. (HI), gives 

/ \ 

+ 00 



j=0 



(7) 



k=0 



In Eq. dV]) the "prime" distinguishes the moments from the central moments. The recursion formula 



that relates the moments to the cumulants 12011 in our case reads 



(„{*)) 



s=l 



[ \ 

J-1 



(8) 



from which it can be shown that the j-th order moment, fi'^iv^^^), is a polynomial of the first j 
cumulants, ^i^' with s = 1, 2, . . . , j. Thus in Eq. ([8]): 



K 



^^2 T 1^*1 ) 



where the coefficients of the different terms are those that occur in the Faa di Bruno's formula. By 
using Eq. ([3]) we can rewrite these terms in the form 

^'^{v^^)) = kn^^'^ + k^iKt"^f 

+iok^K^^'\K^;''^f + k^{K^;''^f 

^^',{v^'^) = . . . , 

in which k is the number of detected photons, each one converted with its own 7^(2 = 1, 2, . . . , fc), 
and the cumulants are those of the probability distribution p^. We observe that each term contains 
a product of cumulants in which the sum of the indices is the order of the moment. 

Let us assume a narrow p-y distribution so that 0. In terms of cumulants this rewrites 

K^2'^ = o[{K^^'^y]. Under this hypothesis also jls/^^ — > 0, that is k^^^^ = o[{K^^'^y], as we can 
write /i3/7^ = {fi3/a^){a^ /^^), where the first factor is the (finite) coefficient of skewness of 
the distribution p^. Actually it can be easily shown that k^j'^^ = o[(k[^''*)-'] relations hold for any 
j > 2, if ^2^'^ = o[(k[^'-')^]. Taking into account that the number of detected photons, k, is a finite 
number, all the monomials in the above expressions of the moments are negligible with respect to 
the last one, so that we can approximate *^^)) = A;'*(k^^'-')* = k^^y^. By substituting in Eq. (|7]) 
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we get 
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VJ 
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Y Y ^rn=kik - {m)Y = YfJ'T 
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(9) 



fc=0 



Note that, as iir{m) never vanishes, even in the case of light in a single-mode Fock state because 
of the non unit quantum efficiency of the detectors, actually Eq. ^ holds for measurements per- 
formed on optical fields with any statistics. 

Dividing both members of Eq. ©by {v) yields 



{v) ^ (m) 
while the exact results for r = 2 and r = 3 would be 



(10) 



(v) 
{v) 
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{m) Y 
/X3(m 



(11) 
(12) 



(m) ' (m) Y Y 
respectively. 

We thus assess that, when detector and processing electronics ensure a sufficiently small ratio 
(j/7, the scaling law in Eq. (fTOl) holds and the simple knowledge of 7 allows reconstructing P^- In 
fact binning the v data of a measurement into bins of width 7 produces a distribution identical to 
Pm- Alternatively we can say that Pm is recovered by dividing the v data by 7 and then rebinning 
the new values into unitary bins. 

How to determine 7 when o" ^ 7 has been already shown [6]. Here we demonstrate that, 
for a detector simply endowed with linear response, we can both determine 7 and decide on the 



negligibility of a with respect to 7. 

At this point, for the first time in this work, we make use of Eq. ([T]). Such a link between the 
Pm and Pn distributions gives (m) = ri{n) and (m^) = rf(n^) + rjil — r]){n). Thus we find |16|] 

^=r,Q + l. (13) 

where Q = [/i2('^) — (^)]/(^) is the Mandel parameter of the light entering the experimental 
apparatus in Fig. [Hand containing {n) photons in the Tm time interval Substituting Eq. (fT3l) 
into Eq. (flOl) and taking into account Eq. (H)) yield 

^^lM= (1 + ^) . (14) 

in which Q/ (n) is independent of r]. On the other hand, {v) depends on r], which can be changed 
by acting on the light delivery optics: r] can be set at any value between the product of r]q times 
the coupling efficiency of the optical delivery system and zero by adding filter into the system that 
delivers the light to the detector. Thus by repeated measurements of the same light at different i], 
we can verify the linear dependence on (v) in Eq. (fT4l) . The experimental /U2(f )/ {v) data plotted 
as a function of (v) should align along a straight line, whose intercept reduces to 7 if cr^/7^ ^ 1. 
Knowing 7 allows proceeding to the rebinning of the v data that leads to the reconstruction of 
xperimental applications to some non-trivial classical states are described in references 



P 



|6L 
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22] 



ni. DISCUSSION 

The assessment that an experimental apparatus has aa/'j ratio sufficiently small for the validity 
of Eq. (flOl) deserves some comments, owing to the difficulty of knowing 7 and a separately. 

We first observe that for any photo-emissive detector a decreases at increasing the strength of 
the electric field experienced by the photoelectrons as soon as they leave the cathode. For a PMT in 
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which the internal gain is provided by multi-dynode cascade amplification, increasing the voltage 
difference between cathode and first dynode produces smaller a values. For PMT's in which the 
electron amplification is provided by other structures {e.g. micro-channel-plate, metal channels, 
etc.), the same effect is obtained by acting on the voltage of the accelerating electrode. For a HPD, 
in which the electrons released by the photocathode are multiplied by a reverse biased avalanche 
diode, cr is lowered by applying greater negative high voltages to the photocathode. 

In the case of PMT's, modifying the voltage partition to change a brings about a change in 
7 that cannot be easily compensated by acting on the overall voltage applied between anode and 
cathode. In the case of HPD's this compensation is feasible by adjusting the avalanche diode 
reverse bias voltage. However in any electronic apparatus that processes the detector output there 
is a step that allows changing 7 by a known factor {e.g. in Fig. [TJ both AMPL gain and ADC 
scale) while keeping (j/7 virtually constant. 

The expression of the intercept in Eq. (fT4l) is such that, upon changing 7 by a known factor 
but not cr/7, a new series of measurements at different rj values would provide a new evaluation 
of the intercept, whose value should scale by the same factor. On the contrary, for constant 7 and 
different (j/7 ratios, the intercept should change differently. Note that a check of the constancy of 
7 is provided by Eq. dH) in which a does not appear. If, by manipulating the voltages supplied to 
PMT/HPD detectors as described to change a, we achieve a situation of constant and minimum 
intercept, we have proved that (T^/7^ ^ 1 in Eq. (fT4l) . We can thus use this limit value of the 
intercept as the correct 7 to rebin the experimental distribution and reconstruct Pm- If the (m) 
value provided by the reconstructed P^ fits Eq. ([H), it means that the detector guarantees a (T^/7^ 
not only much smaller than one but small enough for the validity of Eq. (flOl) . 

We finally note that the above described verifications of the validity of Eq. (flOl) are self- 
consistent in that they do not require measuring to establish the relation between and 7^. 

11 



This is a definite advantage with respect to any potential method for reconstruction based on 
the determination of the moments [18]. 



IV. CONCLUSIONS 



We have shown that for any linear detector we can both measure 7 and determine if cr/7 is 
sufficiently small for taking as reliable the Pm reconstruction achieved by binning the experimental 
V values into bins of width 7. For the method to work it is necessary that the m-range where Pm 
is non-negligible falls within the linearity range of the apparatus, which must be broad enough 



for a satisfactory verification o 
works not only with HPD's [16, 



" Eq. (fT4l) . In forthcoming papers we will show that our method 
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and the Burle 8850 PMT |15il6|], but also with detectors such 



as Si multi-pixel photon detectors [|23ll . and more PMT's endowed with single photon sensitivity. 
Useful detectors might also be solid state detectors such as avalanche photodiodes in the linear 
amplification regime jl^ . By the way, some photon-number resolution is being demonstrated for 
these detectors, in particular if connected to charge-integration readout circuits with sufficiently 
low noise [|25ll . At last, for a thermal detector such as a TES, obtaining a a value sufficiently 
smaller than 7 would be a minimal performance as compared to the excellent photon-number 
resolving power demonstrated by Lita et al. up to 7 detected photons and might allow using a 
less sophisticated apparatus. 

We think that the results described in this paper will broaden the choice of detectors suitable for 
measuring photon statistics. The essential requirement for the detector, beside that of the linearity 
of the response, is the smallness of the ratio cr^/7^, which can be ascertained (see above) without 
measuring p^. The fact that the method applies to measurements in the macroscopic realm may 
turn out to be relevant in all cases in which one cannot attenuate the light to bring the photon 
detection rate down to the regime where photon-counters operate. As examples we mention fields 
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produced by events either rare or unstable and, more importantly, all nonclassical fields, where our 
method risks being the only one applicable to macroscopic fields. 
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